Orbital Dependence of Quasiparticle Lifetimes in SrzRuO^ 
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Using a phenomenological Hamiltonian, we investigate the quasiparticle lifetimes and dispersions 
in the three low energy bands, 7, f3, and a of Sr2RuOi. Couplings in the Hamiltonian are fixed 
so as to produce the mass renormalization as measured in magneto-oscillation experiments. We 
thus find reasonable agreement in all bands between our computed lifetimes and those measured in 
ARPES experiments by Kidd et al. Q and Ingle et al. 0. In comparing computed to measured 
quasiparticle dispersions, we however find good agreement in the a-band alone. 



INTRODUCTION 

The ruthenate SV2-RUO4 has been extensively studied 
for a number of reasons, chief among them its unconven- 
tional superconducting state [2j . Rs electronic structure, 
reviewed in detail by Bergemann et al. , has also been 
of much interest. Three bands belonging to the ti g com- 
plex of 4d Ru orbitals cross the Fermi energy. They are 
divided into two sets. One derived from the d xy orbital 
has a two-dimensional dispersion with little dispersion 
along the c-axis due to the layered structure of the ma- 
terial. The second set comprises the d xz and d yz bands 
which have predominantly one-dimensional dispersion. A 
key feature is the absence of hybridization between these 
sets in a single layer due to the opposite parity under re- 
flection about a R11O2 plane. This contrast in their dis- 
persion has been invoked by Kidd and collaborators [l| to 
explain the strong difference in the energy and temper- 
ature dependence of the quasiparticle lifetimes between 
the two sets observed in recent ARPES (Angle Resolved 
Photoemission Spectroscopy) experiments. In this note 
we report on some simple calculations to examine this 
interpretation. 

The Fermi surface is illustrated in Fig. 1 and consists 
of three sheets. The almost circular 7-sheet derives from 
orbitals with d xy symmetry. The a and (3 sheets derive 
from the approximately linear sections due to the orbitals 
with d xz and d yz symmetry which however hybridize 
weakly with each other where they cross. Kidd and col- 
laborators measured the dispersion and linewidths of 
the quasiparticles on the 7- and the (3- sheets along the 
line F-M and found a clear difference in their behavior as 
a function of both energy and temperature. Ingle et al. 
measured the same quantities for the a-band along 
the F-X direction finding similar linewidths to those ob- 
served by Kidd et al. in the /3-band. At first sight this 
difference between the 7-band and the a-, /3-bands would 
seem to be simply a consequence of their differing disper- 
sion and the lack of hybridization between them. While 
direct scattering of a quasiparticle between these bands 
is forbidden due to their different parities, quasiparticles 
in these band will still interact through the Coulomb in- 
teractions which leads to modifications of their strictly 




FIG. 1: A schematic of the three bands on SriRvOt,. 



one and two-dimensional character. To explore this ef- 
fect we undertake low order calculations of the lifetime 
using phenomenological interaction strengths. These in- 
teraction strengths are chosen so that the mass renor- 
malization observed in magneto-oscillation experiments 
0,0 is reproduced. 

Doing so we are able to account for the broad features 
of the observed linewidths. In particular we obtain a 
good fit to the linewidths in the 7— band as a function 
of binding energy and we find, consistent with experi- 
ment that the linewidths in the (3- and a-bands are far 
smaller. Our computations do however consistently un- 
derestimate to a small degree the linewidths. This points 
to, perhaps, some other contributory mechanism beyond 
electron-electron interactions. We also explored the fi- 
nite temperature behavior of the linewidths. In both the 
7- and (3- bands, our computations match well the ob- 
served behaviour seen in |l|. Unlike the linewidths, the 
agreement between the measured and computed disper- 
sion relations is mixed. For the a-band, good agreement 
is found while for the (3 and 7-bands, the measured mass 
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renormalization is far less than what would be expected 
from magneto-oscillation experiments. We comment on 
this further in the discussion and conclusion section. 



CALCULATIONS OF THE LIFETIME AND 
EFFECTIVE MASS 

The multi-band Hamiltonian that describes a single 
layer is as follows 

H = Ho + Hint ; 

Ho = Y, e ^ k Kk« c ^ 

k(Tfl 

Hmt = ^ V,i U n li(r (i)nua'(i)- (1) 

Here the greek indices, /i, v, are band indices and sum 
over the three bands, a, /?, and 7. The two-dimensional 
7-band, describing orbitals with d xy symmetry, takes the 
form 

e 7 (k) = — .88 cos(k x ) — .88 cos(fcj,) — .56 cos(k x ) cos(k y ) — .5. 

(2) 

The one-dimensional a- and /3-bands arise from weak hy- 
bridization between the d xz and d yz symmetry orbitals 
and appear as 



e xz (k) 



f v z 



— .62 cos(afca;) — .09cos(afcj,) 

+.04 cos(ak x ) cos(aky) 
(k) = — .09 cos(afc x ) — .62 cos(afcj,) 

+.04 cos(ak x ) cos(aky) 

e ± (k) = he"±e"*); 



.24; 
.24; 



*//3 _ P + 



+ V(e-) 2 + -01. 



(3) 



All energies are in eV's and a = 3.86^4° is the lattice 
spacing. The interaction vertices reduce to four terms 
describing intraband scattering between electrons in the 
7-band (V-y-y) and in the a„3-bands (V aa — Vpp) and in- 
terband scattering between 7 and the a, (3-hands(V ia — 
Vyp) and between the a and (3 bands (V^g). Simple esti- 
mates of the Fermi-Thomas screening length give strong 
screening due to the large density of states at the Fermi 
energy so we neglect the wavevector dependence of the in- 
teraction vertices. This leads to four independent param- 
eters which we adjust phenomenologically as discussed 
below. The calculation of the self energy is restricted 
to the lowest orders in the phenomenological interac- 
tion vertices, illustrated in Fig 2. For the lifetime this 
amounts to calculating the decay processes of a single 
quasihole into three quasiparticles with renormalized in- 
teraction vertices. These processes will dominate the de- 
cay rate of low energy quasiholes as multiquasiparticle 
decays will rise more slowly with the energy of the quasi- 
hole measured from the Fermi energy. Our aim is to 
explore the effects of the differing band dispersions on 
the decay rates. 
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TABLE I: Mass renormalization of the electrons in the three 
bands of SrzRuOi as given by magneto-oscillation experi- 
ments 0,Q (dHvA), ARPES measurements QQ, and LDA 
computations . These three mass renormalizations are all 
given in terms of the bare electron mass. 



Our approach to this problem is then in the spirit of 
Landau Fermi liquid theory - we use the necessary mass 
renormalization to fit the coupling constants (akin to the 
Fermi liquid parameters) and then determine other quan- 
tities (the inverse lifetimes) in terms of these same cou- 
plings. In this spirit, the largeness of certain param- 
eters (for example, the dimensionless coupling it 77 = 
KF'yV T y/2'KVFf), should not be taken as problematic. 

The real part of the self energy is calculated employing 
the same set of diagrams. Any non-zero contribution to 
the real part at uj = at the Fermi surface is absorbed 
into a redefinition of the chemical potential. The first 
order diagram, marking the exchange energy, contributes 
to the real part of the self-energy alone. It is the first in 
a series of diagrams but is the only member of the series 
to make a contribution. The remaining diagrams in the 
series simply renormalize the chemical potential and so 
are ignored. 

To compute the strengths of phenomenological inter- 
action parameters we insist that the computed self en- 
ergy reproduce the mass renormalizations expected from 
magneto-oscillation experiments 0, 0] (see top line of Ta- 
ble 1). The mass renormalizations induced by the inter- 
action terms in Eq. 2.1 are given by 



m R 



1 - 9 w ReS(w,A: = K F ) 



mLDA 1 + v L p A dkReY<(w, k = K F ) 



(4) 



where tulda is the mass renormalization induced by the 
LDA band structure (also listed in Table I). The LDA val- 
ues given in Table I are averages around the bands' Fermi 
surfaces weighted by the local density of states (appro- 
priate as the magneto-oscillations masses are themselves 
weighted averages). The interaction strengths are then 
determined by insisting that trr equal m^/mLDA (as 
listed in Table 1) at a single point along the Fermi sur- 
face of each band. We do not attempt to average wir 
itself around the Fermi surfaces (requiring a computa- 
tion of the self-energy everywhere) as too computation- 
ally costly. 

In the presence of contact interactions, the T = self- 
energy for a given band fi at first order consists of a term 
independent of both momentum and energy and so may 
be ignored (at finite T, this term however does do more 
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FIG. 2: Diagrams contributing to the self energy at first and second order. Here the greek indices a, j3 refer to generic bands. 
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TABLE II: Strengths of the interaction couplings (all values 
in a 2 ■ eV where a is the lattice spacing. 



than renormalize the Fermi energy) . At second order the 
contribution takes the form 



= E^/ rf « j dx 



w + i5 — e^(k — q) — x 
x(e(-x)-6(e,(k-q))), (5) 

where \v is the susceptibility of electrons in band v. We 
numerically compute the derivatives of the real part of 
the self energy necessary to computing run in order to 
accurately capture the effect of the dispersion in each of 
the bands (without recourse to approximating the bands 
through a quadratic dispersion relation). 

The values of the interaction parameters so determined 
are given in Table II. The first number marks the value 
used in later computations of the dispersion and lifetime. 
The ranges in brackets mark the range of parameters pro- 
ducing the desired mass renormalization, m^. However 
other choices within the given ranges yield results for the 
lifetimes that differ by no more than 10%. The largest 
coupling by far is V~ n , dictated by the need to produce 
a mass renormalization of mji ~ 7. 

The real part of the self-energy will have terms beyond 
linear order in uj that however are not necessarily negli- 
gible a priori. In three dimensions the real part of the 



self energy would be expected to take the form 

ReE(w) = aui + bu; 3 H . (6) 



By dimensional analysis b would of the form, 

J F 



dimensionless constant x E F 2 , and so for uj small, the 



cubic term can be ignored. But in two dimensions, the 
self-energy develops non-analyticities and will appear as 



ReS(w) = aui + buj 2 sgn(uj) + 



(7) 



It thus cannot be immediately neglected. 

We can estimate this non-analytic term more precisely. 
Using jjj], we know the non-analytic contribution to the 
self-energy at second order at arbitrary uj and k is given 
for bands with quadratic dispersions [lOj by 



ReS^° n ~ anal -(w,fc) = 



V 2 m 2 



(uj 2 + -(u - e M (fc)) 2 )sign(w - e p (*)) 

+ {lj 2 - -(u - e M (fc)) 2 )sign(u; + e M (fc)) 
+cj 2 sign(w) + w 2 sign( — — ) 



-(w + e M (fc)) 2 sign(- 



— anal. 



V 2 v m^ml 
2>2^Kp^K Fv 

\2\ 



W + 4^ - £ mO)) )sign(w - e M 



(*)) 



+ (uj 2 - -(u - e M (/c)) 2 )sign(w + e M (fc)) 



(8) 



where E Ffi is the effective bandwidth of band /i, and 
rrifi = K F /2E Ffl with K Ffl the band's corresponding 
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FIG. 3: Plots of the inverse lifetimes in the three bands, 7, /3 and a. 



Fermi wavevector. £^° n ~ anal -(ti>, k) has two contribu- 
tions: one arising from the non-analyticity in the elec- 
tron polarizability, x^O^i*?)) a t Q near and one from 
another non-analyticity in q near 2Kp^ 0- I n contrast, 
the interband contribution to the self energy, E M „, only 
sees a non-analytic contribution from ^(a;,^ ~ 0). To 
estimate the contribution these terms make to the self- 
energy we need to estimate the appropriate values for the 
effective (quadratic) masses appearing in the above. To 
do so we write down the corresponding expressions for 
the imaginary part of the self energy: 
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FIG. 4: Plots of the dispersion in the three bands, 7, f3 and a. 
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TABLE III: The effective bandwidths for each of the three 
bands. 



Given that u 77 ~ 5 and Ep-y ~ 0.5eV, for very small u> 
we would then expect higher order contributions to lead 
to a reduction in the inverse lifetime while for omega 
larger than .02eV, higher order contributions leave the 
answer effectively unchanged. 



To determine Ep^, we fit a numerical evaluation of the 
imaginary part of self-energy via Eqn. 2.5 to the above 
expression. (The non-analytic portion of ImS(w, k) dom- 
inates because of the presence of the logs.) Having so ex- 
tracted effective Fermi energies, Ep^ (see Table III), we 
then are able to evaluate the non-analytic real part of the 
self-energy. We find that it makes a significant contribu- 
tion only for ui > .04eV and then only for Re£ 77 - the 
dispersion relations for the a and /3-bands are insensitive 
to this correction. 

While the spirit of our computation is Landau Fermi 
liquid theory, we can nonetheless address the question of 
what contribution higher order diagrams (beyond those 
in Figure 2) make to the inverse lifetime. This is a dif- 
ficult question of course, but for at least an infinite sub- 
set of diagrams studied by Chubukov et al. |{J, we can 
give a definitive answer: not much. We focus on E 77 
(i.e. the contribution to the self energy of the 7-band 
due to the polarizability of 7-band electrons) as V 71 is 
the largest coupling. Chubukov et al. have studied 
the behaviour of higher order forward scattering terms 
in the V^ 7 -perturbative series. They found that the rel- 
evant expansion parameter is s — u 2 ^u>/(uj — e 7 (fc)). In 
terms of our lifetime computations, we are far off the 
mass shell because of the large mass renormalizations in- 
volved, i.e. (uj — e-f(k))/u> ~ 6, s will be smaller than 
might first appear, i.e. s ~ 4. While the radius of con- 
vergence of the series is unclear, if the series needs to be 
resummed, the effects are not particularly drastic. The 
terms in ImE(w, k) (see Eqn. 2.8) behaving as uj 2 ln(cj) 
are modified to 

u? ln(-^) - c 2 (l ln(-£_) - h \n(u 2 )\) (10) 



COMPARISON TO ARPES RESULTS 

The ARPES studies of Kidd and collaborators Q re- 
ported inverse lifetime results for a hole in the 7 and /3 
bands measured along the T — M direction as a function 
of energy and of temperature. There is a pronounced dif- 
ference between the bands with the inverse lifetime for 
holes in the 7-band rising much more rapidly with energy 
and temperature than in the /3-band. The results for the 
/3-band are mimicked by the linewidths observed in the 
a-band in 0. 

In the top panel of Figure 3, we plot our computed 
inverse lifetime as a function of binding energy for the 
7 and /3-bands along T — M. To our computed result 
of the 7-band, we have added a static impurity contri- 
bution of ImEimp = .055eV. This is chosen so that the 
experimental and theoretical values of the inverse lifetime 
coincide at uj = 0. To arrive at the experimental values 
of the inverse lifetime, we have taken the MDC widths 
as measured by Kidd et al. and multiplied them by the 
corresponding LDA value of the band velocity, vf-lda- 
For the 7-band this velocity changes rapidly near the 
Fermi surface and so we have correspondingly employed 
a fc-dependent vf-lda- This yields a slightly larger ex- 
perimental inverse lifetime than reported in While 
our computed values undershoot the measured values of 
the inverse lifetime, the agreement is reasonable for the 7- 
band and would be improved significantly if we treated 
ImEimp as a completely free parameter. In particular, 
we note that the uncertainties in the measured values do 
not reflect the uncertainty in the correct value of the bare 
LDA velocity to employ in converting an MDC width to 
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FIG. 5: Plots of the temperature dependence in the 7- and /3-bands. 



an inverse lifetime. 

For the /3-band the computed and measured values of 
the inverse lifetime broadly match. In particular, the pre- 
dicted and measured linewidths are much smaller than in 
the 7-band. Again we have added an impurity contribu- 
tion to the computed inverse lifetime values so that the 
measured and computed values agree at zero binding en- 
ergy. We do note however that the computed inverse 
lifetime is smaller than that measured. Over a range of 
binding energies of .06eV, the measured inverse lifetime 
increases by ~ .05eV whereas the computed value in- 
creases only by half that. This disagreement may be in 
part explained by the presence at larger energies of cor- 
respondingly large error bars. However, unlike for the 
7— band, the band velocity in the /3-band is not a strong 
function of wavevector. 

In the bottom two panels of Figure 3, we plot our com- 
puted lifetime for the a-band along the M — X direction 
vs that measured by Ingle et al. Q. Ingle et al. use 
samples with freshly cleaved surfaces where the a-band 
can be resolved without employing surface aging. In the 
left panel, we plot our computations, again with a cor- 
recting Im£j mp , against the values measured by Ingle et 
al.. In the right panel we plot our computations (with 
no impurity correction) against the measured electron- 
electron contribution to the self-energy. Ingle et al. de- 
duce this contribution by subtracting out an estimate of 
both the phonon contribution and the impurity contri- 
bution. Again we see that the computed value of the 
inverse lifetime in the a-band is much less than that of 
the 7-band. However like the /3-band, this value is less 
than that measured. If we look at only the value imputed 
to the electron-electron contribution, we see a rise in the 
inverse lifetime of ,05eV over a .06eV range of binding 
energies whereas our computed value rises by a little less 
than .02eV. 

In Figure 4 we compare the corresponding computed 
and measured dispersion relations for the three bands. 
The computed dispersions are essentially straight lines 



whose slope is determined by the value of the dHvA 
mass renormalizations (as this is how we determine the 
coupling strength). We have also plotted the LDA dis- 
persion relations. For the 7- and 0- bands, there is 
strong disagreement between measured and computed 
values. Alternatively, the measured dispersion for the 
[3- and 7- bands show a mass renormalization far smaller 
than seen in magneto-oscillation experiments. For the 
/3-band, the measured mass renormalization equals its 
LDA prediction. In the a-band, in contrast, the disper- 
sion closely matches the computed prediction, that is, the 
mass renormalization in the a-band measured by ARPES 
closely matches that measured by dHvA. 

Finally in Figure 5 we compare the computed tempera- 
ture dependencies of the inverse lifetimes at zero binding 
energy in the 7- and /3-bands to that measured by Kidd et 
al. [l|. The /3-band results are well-matched. Both com- 
puted and measured values show only a weak dependence 
on temperature over a 100K range. Our computations, 
nonetheless are consistent with a quadratic Fermi-liquid 
like dependency on the temperature. 

The 7-band results also agree well with that reported 
in [l| . We roughly expect the temperature dependency to 
behave as AtT 2 \og(Ep /T) where At is some constant. 
Similarly we expect the dependency of t _1 on binding 
energy to be (at leading order, modulo the complications 
of Eqn. 9) A u u> 2 \og(Ep /oj). From fits of the two com- 
puted curves we find that At/A w ~ 10. The value of this 
ratio compares well with the expected value, 7r 2 , from a 
single quadratic band @. 



DISCUSSION AND CONCLUSIONS 

A primary conclusion from our analysis to the inverse 
lifetimes and effective masses is that there is a marked 
difference between the two sets of bands, 7, and a and (3 
in terms of effective coupling strengths. This difference 
is not a simple consequence of the band dispersion and 
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its effect on the density of states in intraband scattering 
processes, but rather it is due to a stronger effective in- 
traband interaction in the 7-band at low energies. This 
conclusion is rather surprising since all three bands de- 
rive from orbitals belonging to the same ti g manifold and 
standard multiband Hubbard models assign equal onsite 
interactions to all three orbitals 0, • O ur results point 
to substantially different renormalizations of the low en- 
ergy effective intraband interactions in the two sets of 
bands. 

The different standing of the band 7 from the bands 
a and (3 may arise from the latter two band's highly one 
dimensional character. This character raises the ques- 
tion of whether some form of Luttinger liquid behavior is 
occurring. The fact that the transport properties show 
clear three-dimensional Landau- Fermi behavior shows 
this behavior dominates at low energy. Nonetheless a 
crossover to Luttinger liquid behavior could appear at 
higher energies. However the ARPES data reported by 
both Kidd et al. Q and Ingle et al. 0] show no sign 
of such a crossover up to energies of 60 meV. Both the 
P and a bands appear to have inverse lifetimes that de- 
pend quadratically on the binding energy. This perhaps 
is due to the hybridization and interband scattering be- 
tween the three bands which act to suppress a strictly 
one-dimensional character in the a, (3- bands. 

Our computed inverse lifetimes broadly match the 
scale of the observed linewidths in all three bands. 
Nonetheless they also consistently underestimate the cor- 
responding measured values, more pronouncedly in the 
a- and (3- bands than in the the 7-band. This might sug- 
gest that some additional mechanism is making a con- 
tribution to the self energy, at least in the a- and (3- 
bands. This might point to a possible role played by 
some bosonic mechanism such as phonons. However if 
the shortfall is to be explained by phonons, phonons need 
to make a greater contribution to the self energy than 
posited in Ingle et al. Q{. As we can see from the bot- 
tom right panel of Figure 2, the portion of the self-energy 
in the a-band measured by 0] due to electron-electron in- 
teractions is double that computed in our effective model. 

In the 7- and /3-bands we see that ARPES predicts 
a mass renormalization far smaller than that found in 
magneto-oscillation experiments. This might suggest 
that there is an effective scale in the problem below 
the sensitivity of typically ARPES measurements (i.e. 
•C lmeV). However we also found that the ARPES dis- 
persion in the a-band closely matches that predicted by 
dHvA measurements. This might suggest that the sur- 
face aging performed by Q to distinguish the bulk 7- 
and (3— bands from the surface counterparts changes the 
mass renormalization in some unexpected fashion. How- 
ever our match to the a-band dispersion measured in |2j 
is not without difficultly. Our match to their dispersion 
is predicated solely on electron-electron interactions. If 
we were to ascribe a role to phonons, it would mean we 



have overestimated the coupling strengths, V^ v . This in 
turn would mean we have overestimated the inverse life- 
times. But our inverse lifetimes are already smaller than 
the measured values. It is unclear whether phonons could 
then self-consistently make up the difference. Nor is it 
unproblematic to have phonons be the dominant inter- 
action in a material where it is believed electron-electron 
interactions are responsible for its superconductivity. 

Of course, the discrepancies found in regards to the 
real part of self energy might simply point to a need to 
go beyond our use of low order diagrams based on an ef- 
fective Hamiltonian. It would be interesting to attempt 
a more sophisticated treatment of the set of interacting 
bands in Sr2RuO^. One approach may be to adopt a 
functional RG approach [HJ . Extensive work of this sort 
has already been done on one-band two dimensional sys- 
tems It should be possible to extend this work to a 
multi-band system. 
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